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Within dynamical mean-field theory, we study the metal-insulator 
transition of a twofold orbitally degenerate Hubbard model as a function 
of a splitting A of the degeneracy. The phase diagram in the — A plane 
exhibits two-band and one-band metals, as well as the Mott insulator. 
The correlated two-band metal is easily driven to the insulator state by 
a strikingly weak splitting A ^ of the order of the Kondo-peak width 
zW , where z ^ 1 is the metal quasiparticle weight. The possible rele- 
vance of this result to the insulator-metal transition in the orthorhombic 
expanded fuUeride NH3K3C60 is discussed. 



I. INTRODUCTION 

Strong electron correlations in multi-band, orbitally degenerate systems represent 
an important current theoretical challenge. A lively experimental playground for that 
is provided by electron- and hole-doped fuUerene systems, which exhibit a variety of 
behavior, including unconventional metals like cubic CsCgo 0] , superconductors of the 
A3C60 family (A= K,Rb,Cs) |||,^, superconducting FET devices and insulators, 
presumably of the Mott- Jahn- Teller type 1^,0] • Among them we can place A4C60 
g, NaaCeo § and the class of compounds NH3K3C60, NHgKaRbCgo, NHgKRbaCgo, 
as well as NHaRbaCgo pO|-p^. In these ammoniated compounds, insertion of the 
electronically inert NH3 molecules acts to expand the Ceo lattice, turning the cu- 
bic, metallic and superconducting state of K3C60 into an orthorhombic narrow-gap 
antiferromagnetic insulator ||lo| , [ll| , |l3|] . In this transition, the increase in volume per 
Ceo molecule relative to K3C60, with its probable slight decrease of electron effective 
bandwidth W , is believed to play an important role, as confirmed by the pressure- 
induced reversion to a fully metallic and superconducting - while still orthorhombic 
- phase |l|Jl|. 

Nonetheless, a possible role of the crystal-structure anisotropy in this Mott transi- 
tion - a transition generally believed to occur more readily for reduced orbital degen- 
eracy 1 16 1^, cannot be excluded. It has even been proposed that precisely the split- 



ting of degeneracy [II8I induced by the orthorhombic distortion could be the crucial 
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ingredient driving the transition from the metal and superconductor to the insulator 
PH l . Earlier work showed the Mott transition in d-orbitally degenerate lattice models 
to take place at larger values of U /W for larger degeneracy, roughly proportionally 
to Vd [ p^lpTf l . A strong enough "crystal-field" splitting of the threefold degenerate 
tiu molecular orbital of Ceo caused by orthorhombic anisotropy could remove one or 
several bands away from the Fermi level, effectively reducing the orbital degeneracy 
d, thus shifting the Mott transition to a smaller critical value Uc/W. Exploring this 
concept, the basic question is how strong this splitting must be to promote such an 
effective reduction of degeneracy. In a non-interacting system, that reduction would 
clearly require a splitting magnitude similar to or larger than the full electron band- 
width W. In a strongly-interacting system, it is important to understand whether 
influencing the metal-insulator transition will again require anisotropic splittings as 
large as the bandwidth, or else if some smaller energy scale will emerge in its place. 

In this paper we exploit dynamical mean field theory (DMFT) to answer this 
question, studying the effects of a band splitting on the Mott transition of an or- 
bitally degenerate, strongly-correlated metal. For that purpose we adopt the simplest 
idealization, namely a, d — 2 orbitally degenerate Hubbard model with bandwidth W 
and on-site repulsion U , and mimic the noncubic anisotropy by a diagonal splitting 
A between the two bands. In this model, we study the T = phase diagram in the 
{U, A) plane, and examine in particular how the Mott transition is influenced by the 
splitting. As the aim of the calculation is understanding a simple, and possibly fairly 
general physical mechanism, we do not attempt a realistic calculation of the kind that 
was accomplished recently for Ni and Fe [ pO[ . 

The main result found is that the effect is remarkably strong at large U, where the 
anisotropy-induced transition between the correlated multi-band metal and the Mott 
insulator takes place very readily, in fact for A ^ W. We interpret that as follows. 
The metallic band, originally of width W , is renormalized in the strongly correlated 
metal into a narrower "Kondo-like" quasiparticle peak of width zW , with 1 3> z — > 
as U Uc- Close enough to the isotropic and orbitally degenerate Mott transition 
at U ^ Uc, the correlated two-band metal can be driven effectively one band by a 
splitting as small as the Kondo energy scale A cx zW. The one-band metal, in turn, 
has a lower Uc, so that the resulting state is an insulator. 

With this model result in hand, we move on to consider NII3K3C60, where we 
carry out a DFT-LDA (density functional theory in the local density approximation) 
electronic-structure calculation to extract an order of magnitude for the effective 
orthorhombic splitting parameter A/W, found to be of order 0.25. Despite the obvious 
differences between the two-band Hubbard model and the true three band system, 
we suggest that the orthorhombic splitting will similarly drive the initially cubic 
fuUeride with a small quasiparticle residue z from metal to insulator. In turn, in the 
light of our simple-model phase diagram, the return to metallicity of the ammoniated 
orthorhombic fuUeride under pressure leads to an interesting question regarding the 
effective number of bands in that metal. It is finally concluded that anisotropy is 
likely a relevant element in promoting the Mott state in this material. 

This paper is organized as follows: Sec. |l| introduces the model and deals with its 
simple limiting cases. The calculation and results within the DMFT are described is 
Sec. HI, The ti„ band structure of NH3K3C60 and its splitting due to anisotropy is 
described in Sec. IV. The discussion and conclusions are finally presented in Sec. ^ 
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II. A SIMPLIFIED TWO-BAND HUBBARD MODEL 



We work with the simplest two-band Hubbard model, assuming purely diagonal 
hoppings between orbitals on different sites, the anisotropic symmetry lowering all 
embodied in a diagonal on-site splitting term. We write the Hamiltonian as: 

H ^-t E ("^LTCjacT + H.C. ) + E E '^^a + ^ E ~ ^) ' 

(ij-),cr a=l,2 i Q=l,2 i 

where c]^^ creates a spin-cr electron in orbital a at site i, nia — X]<t '^iao^iaa, n-i = 
X^a is the total number of electrons at site i, and denotes nearest-neighbor 
sites. The first term is a standard tight-binding hopping, simplified to ignore non- 
diagonal hopping, and also merohedral disorder, both of them present in the real 
systems. The second term introduces the anisotropic splitting A as a shift of the on- 
site orbital energy £2 = — ei — A/2. The last term is the Hubbard on-site interaction, 
which we take in its simplest form, omitting for example any additional intra-site 
Hund's rule interaction terms. We also ignore the on-site electron-vibron couplings 
and the ensuing Jahn- Teller effects, even if quite important for other aspects |^,^,^. 

In the U = limit, the splitting A > simply operates a rigid shift of band 2 
upwards with respect to band 1, promoting electron transfer from the upper to the 
lower band. Above a critical value A = Ac, and for a total electron filling n < 2, 
the upper band is emptied up. For example, with two symmetric bands of width W 
and n = 1 electrons per site, the upper band is emptied above Ac/W = 0.5. At zero 
temperature (where we shall restrict ourselves in this paper) and zero interaction, 
[/ = 0, the transition between the "two-band metal" and the "one-band metal" is 
continuous. Because the topology of the Fermi surface changes, this transition is 
accompanied by a weak singularity of the total energy first described by Lifshitz [ p3[ . 
When the electron-electron interaction U is turned on, one expects the emptying of 
the upper band to take place at smaller values of Ac, owing to the effective band 
narrowing. Perturbatively in U one can show that 

A,{U) ^ A,{0) ~ J U + 0{U^) , (2) 



where the value of the coefficient 7 > depends on details of the bands. In addition, 
electron-electron interactions might modify the nature of the metal-metal transition 
singularity relative to the non- interacting case , a point which we will not further 
address here. 

We restrict to a filling of one electron per site, n — I, discarding the trivial case 
n = 2 where the large-A state is a band insulator, as well as noninteger n, where 
insulators are not possible. The A = and A ^ cx) limits reduce then, respectively, to 
the (quarter-filled) two-band and the (half-filled) single-band Hubbard models, both 
possessing a metal- insulator transition as a function of U [p5|-p8t. In particular, we are 
not interested in the weak-coupling antiferromagnetic instability, of the ensuing large- 
A half-filled band, associated to nesting in specific hopping structures: we assume 
a genuine Mott transition for the half-filled single-band model to occur at a finite 

jjl-band > 

The limit of strong interaction U ^ W is insulating for any value of A. This 
limit can be studied by mapping the model (^ onto a spin and orbital exchange 
Hamiltonian which reads |^9[] : 

i?exch = J 51 (S»-Sj+T'-Tj+4S,-S, T,-T,)-A^7;\ (3) 

<ij> i 
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where the pseudospin-1/2 operators T^- = 1/2 J^afii' ^jau'^aP'^jPf represent the orbital 
degrees of freedom, and J = 2t^ /U . For A = this model has been studied both 
in one |30[| and two dimensions with suggestions that interesting spin- liquid 

physics could be realized. For our purposes, it suffices to note that the model has 
no ferro-orbital instability, and has therefore a finite q — orbital susceptibility. As 
a consequence it takes a nonzero value of A to fully orbitally polarize the ground 
state. Due to the absence of cross-band terms in the kinetic energy, complete orbital 
polarization occurs at a finite Aj, oc J = 2t^ /U . For A > A^, that ground state is 
represented by a one band Mott insulator plus a totally empty split-off band. 

Notice, finally, that for the specific case we are interested (n = 1 electron in a two- 
fold degenerate band) no orbital ordering is present at weak coupling; the possibility 
of antiferro-orbital ordering within the Mott insulating phase, on the other hand, 
depends crucially on the details of the hoppings and lattice coordination, and it is, in 
any case, beyond the scope of any single-site DMFT study. 



III. PHASE DIAGRAM OF THE SIMPLE MODEL 

We have now all the elements to illustrate the general layout, sketched in Fig. |^, 
of the zero-temperature phase diagram of model in Eq. (|l]) as a function of ([/, A) 
for n = 1 electron per site. The metal-insulator transition for n=l, d = 2 and 
A = 0, previously studied in DMFT by quantum Monte Carlo was found at 
jj2-band ^ 1,5 W at the relatively low temperature T — W/32. The same transition 
with n = 1 in the one-band model (d = 1, or d = 2 at A = oo) has been similarly 
calculated at [/i-band ^ ^.3 also at T = W/32 |2^. Finite temperature of course 
must affect somewhat these numerical values, but a common feature with our zero- 
temperature DMFT calculation described below, is that ;72-band > [/i-band_ rj.^^ 
line in Fig. |l|, representing Uc as a function of A and separating metals from Mott 
insulators, indicates this trend. The DM' line separates the fully orbitally polarized 
Mott insulator (at the right) from the two-band insulator, roughly as Uc ^ A~^ for 
small A. Similarly, the CM line separates the fully orbitally polarized metal from 
the two-band metal: it starts from point C with a linear slope —7^^, according to 
Eq. (H). In the region of full orbital polarization the value of A is irrelevant, and 
this is the reason why the Mott transition line MB is horizontal. As A increases, for 
U < C/^^'^^"'* the upper-band emptying transition takes the two-band metal across 
the CM line over to a one-band metal, while for U > [/i-^and leads across the AM 
line to a Mott insulating state. The AM line starts off linearly at small A. Basically, 
this occurs since in close proximity of the Mott transition the relevant energy scale, 
to be matched by the orbital splitting A, is no longer the bare bandwidth, but rather 
the width of the Kondo-like quasiparticle peak of the correlated metal. That in turn 
is of the order of zW, which is smaller than W and proportional to the quasiparticle 
residue z, linearly vanishing with ([/2-band _ jj^ ^j^g Mott transition [|9|. As 
the DMFT calculations will show, the effective emptying transition occurs when A 
increases to reach Ac{U) oc zW oc (^/s-band _ jj-j^ 

The purpose of the DMFT calculations is to investigate this phase diagram in 
greater quantitative detail, in particular in the region of intermediate U ^ W and 
A < M^. We have implemented DMFT in the exact diagonalization flavor [ p^j27| 
for the degenerate Hubbard model on the Bethe lattice, restricted to paramagnetic 
states . The advantage of the exact diagonalization method over the Monte Carlo 
approach in the solution of the impurity problem rests primarily in its being less 
computationally intensive - useful when studying the phase diagram of a model as a 
function of more than one parameter - and in its ability to address zero-temperature 
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properties The results of the DMFT calculations are summarized in Fig. |[ 

The values we obtain for J/i-b^nd ^ 1.35;^ and t/^-i^'^'^d ~ 1.8 VF, associated 
to the disappearance of the metallic state for one and two bands respectively, have 
been calculated with the same number of discretized bath levels per band rif, — 4. 
The associated error bar is estimated by repeating the calculation with rih ~ 5. 
These values are in substantial agreement with corresponding values obtained by other 
methods 1^5^281]. The other points in the phase diagram are obtained by following 
the stability of the two-band metal for a given value of U and increasing A, marking 
the emptying transition to the one-band metal or to the insulator. 

A deficiency of this single-site DMFT calculation - which is restricted as usual to 
paramagnetic states only - is the absence of a two-band insulating state for A > 0. 
In fact, the suppression of the antiferro orbital fluctuations embodied in the exchange 
model (H) produces a fictitious infinite uniform orbital susceptibility which leads to full 
orbital polarization as soon as A is turned on Despite this limitation, the results 
of the present calculations are suggestive, revealing the announced sharp reduction of 
the metal-insulator C/c(A) for small but finite A. Indeed, our calculations give a Ac{U) 
which is roughly proportional to the quasiparticle residue z{U) of the undistorted 
(A = 0) correlated metal at U < C/^-band^ 

^ « f3<U) , (4) 

with a proportionality constant /3 ~ 0.3. Since z{U) vanishes as U ~> [/2-band^ most 
likely linearly in ([/2-band _ jj^ jj^^ even a very small A/W can be enough to cause 
a metal-insulator transition in the strongly correlated metal. For example, following 
the bold arrow at U = 1.5 W in Fig. ^ a A value as small as 0.08 W is sufficient to 
cross the transition line from the metal to insulator. 

Further illustrating that, in Fig. |^ we show the behavior of the spectral density 

Ga{io) being the one-particle Green's function of band a, on both sides of the metal- 
insulator transition. The asymmetry of the upper-band spectral density ^2(0;) (solid 
lines) is very pronounced, as this band is nearly {A/W — 0.07) or completely {A/W — 
0.08) empty. As soon as the Kondo-like peaks of the two bands differ enough in energy 
to induce the emptying of band 2, the lower-band spectral density Ai{lu) takes the 
symmetric shape, characteristic of the half-filled one-band Hubbard model. Here the 
Kondo peak disappears completely, as this value oi U > {/i-band p^^g ^^j^g Hubbard 
model of band 1 well inside the insulating regime. 

IV. ANISOTROPIC BAND SPLITTING IN ORTHORHOMBIC NH3K3C60 

The above model calculations show that a small splitting A ex. zW of the orbitally 
degenerate band can drive the metal-insulator transition. We now wish to explore 
the implications that this result - if assumed to be more general than the simple 
model where it was derived - can have on the metal-insulator transition which takes 
place between isoelectronic K3C60 and NH3K3C60 (the former cubic and the latter 
orthorhombic) , and on the insulator-metal transition of NH3K3C60 itself under pres- 
sure. 

To obtain a quantitative estimate of the actual strength of anisotropy in the am- 
moniated fuUeride, we carried out comparative DFT-LDA electronic structure cal- 
culations for K3C60 and a simplified but meaningful model of NH3K3C60 using the 
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PWSCF package based on plane waves and ultrasoft pseudopotentials Q. We took a 
lattice constant a = 14.2 A for fee K3C60, and a = 14.89 A for NH3K3C60, the latter 
with a centered tetragonal unit cell with c/a = 0.91, neglecting the exceedingly small 
difference between a and b. In both calculations a single-Cgo unit cell was used, as 
pictured in Fig. ^ In these model structures and calculations, we ignored merohedral 
disorder, as well as the rich antiferro-rotational structure recently discovered in actual 
NH3K3C60 ||. 

We used 27 Ry and 160 Ry plane wave energy cut-offs for the wavefunctions and 
the charge density, respectively. The Fermi surface was smeared with parameter of 
0.3 eV and a 2 x 2 x 2 Monkhorst and Pack mesh The density of states (DOS) 
was calculated as a sum of Gaussians of width a = 0.03 eV sampling fc-space with a 
uniform 6x6x6 mesh. 

Fig. ^ shows the ti„ bands of NII3K3C60 (a), compared to those of K3C60 (c); the 
corresponding density of states are compared in panels (b) and (d). We find that 
the presence of the inserted NH3 molecules modifies only weakly the essentially pure 
Ceo tiu conduction band, as expected. The difference in the two band structures 
are almost entirely due to the different lattice structures of the molecular centers. 
Small changes appear in the details of the bands throughout the Brillouin Zone, but 
both compounds are predicted in LDA to be three-band metals, their overall tiu 
density of states not dissimilar in the ammonia-intercalated and pristine compound. 
In particular, as apparent from Fig. |^(b), the bandwidths of NII3K3C60 and K3C60 are 
close, W ~ 0.6 eV. The main difference in the two band structures is a splitting of the 
threefold degenerate ii„ band of K3C60 at the F point of NHsKsCeo- This splitting is 
precisely a measure of the strength of the non-cubic crystalline environment seen by 
the tiu orbital on each fuUerene molecule in the orthorhombic structure of NIl3K3C6o- 
Its magnitude is roughly a quarter (0.2 - 0.3) of the total bandwidth; this represents 
the main result of the DFT calculation. If this material were an uncorrelated metal, 
this splitting would have no major consequences, just a change of shape of the Fermi 
surface. The consequences can be much more important due to strong correlations, 
as discussed below. 



V. DISCUSSION AND CONCLUSIONS 

The F-point tiu band splitting - see Fig. ^(a) - can be taken as a crude estimate 
of the value of A ~ 0.15 eV in NH3K3C60, corresponding to a dimensionless ratio 
A/W ~ 0.25. This value is of course too small to determine a complete emptying 
of one of the bare, uncorrelated bands, and indeed all three bands of NH3K3C60 
still cross the Fermi energy Ep. However, correlations play an important role in 
fuUerides, where the value of U is believed to lie between 1 and 1.6 eV thus 
substantially larger than the bandwidth W. As the non-ammoniated fuUerides A3 Ceo 
are still generally in the metallic and superconducting phase, it is not unreasonable 
to assume that their value of U/W p7|-^ is in all likelihood only marginally smaller 
than the [/3-band ^j^g Mott transition in the half-filled d = 3-band cubic system. 
Indirect evidence that ordinary, uncorrelated electronic bands do not exist, and that 
the electron spectral function at E-p can more reasonably be assimilated to a narrow, 
dispersionless Kondo-like resonance was provided for K3Ceo by recent photoemission 
data |§. 

The analysis of the metal-insulator transition in our simplified two-band model sug- 
gests a mechanism for an anisotropy-induced Mott transition in a correlated multi- 
band fermionic system. Clearly, many ingredients which we have left out of our 
calculation are expected to be quantitatively important for a realistic description of 
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the tiu doped fuUerides. However, we suggest that the mechanism which we have 
described in Sees. || and [II, whereby a smaU orthorhombic spHtting of the bands is 
capable, alone, of turning a strongly correlated metal into a Mott insulator, could be 
a qualitatively important ingredient in the description of these systems. Specifically, 
we suggest that, in the appropriate ([/, A) phase diagram (sketched in Fig. ||), the am- 
moniation of K3C60 into orthorhombic NH3K3C60 could correspond to a displacement 
jjlj l similar to that indicated by an arrow in Fig. ^. Inspection of Fig. || shows that, 
if U is close enough to Uc, even a value of A substantially smaller than that we have 
estimated for NH3K3C60 could suffice to drive that metal insulator transition. Very 
recent conductivity measurements |^ on the class of compounds NH3K3_3.Rb3;C6o 
suggest indeed the orthorhombic distortion as the main ingredient driving the Mott 
transition in these systems. 

We should stress that our picture does not exclude a role, in the metal-insulator 
transition, for the accompanying ammoniation-driven expansion of the lattice. That 
role is experimentally proven by the pressure-driven insulator-metal transition, where 
NH3K3C60 is transformed into a metal and a superconductor, despite the permanence 
of the orthorhombic structure |p^ . However a new, interesting, question raised in 
Fig. ^ is what kind of metal does one recover in NH3K3C60 under pressure, whether 
a three-band metal like the cubic K3C60, or a new effective two-band metal. The 
former might in effect require a smaller deformation. It would be very interesting 
in this regard if the Fermi-surface topology of NH3K3C60 could be studied under 
pressure. 

Another interesting - even if perhaps not practically straightforward - test of this 
overall picture could be obtained by applying uniaxial stress to cubic superconducting 
fullerides of the A3 Ceo family. Contrary to standard tendency of hydrostatic pressure 
toward metallization and lower superconducting Tc, the orbital splitting associated 
with the appropriate uniaxial strain could drive some of these compounds perhaps 
first toward higher Tc, and eventually Mott insulating. A similar suggestion has been 
recently put forward, in the context of field-induced superconductivity MM, by Koch 
||4^ , although on rather different grounds. Also, recent calculations ||4^ , |45| of the size 
and the effects of the charge-inducing electric field on the 2-dimensional Cgo layer, in 
the field-induced geometry, suggest that this field induces splittings of the molecular 
levels which might in turn affect the electrons and holes bands in a qualitatitevely 
similar way to the one we propose here for the orthorhombic field in 3-dimcnsional 
fullerides. 

In summary, we have investigated the role of a symmetry-breaking orbital splitting 
on the transition to a Mott insulating state in an initially orbitally degenerate strongly 
correlated multi-band metal. The calculations, carried out in a highly idealized d — 2 
model, have unveiled the Kondo quasiparticle energy scale as that which the non- 
cubic splitting magnitude must reach in order to induce the Mott transition even 
at constant bandwidth and effective Coulomb repulsion. The likely relevance of this 
main qualitative finding to the fullerides K3C60 and NH3K3C60 has been highlighted. 



ACKNOWLEDGMENTS 



We are indebted to O. Gunnarsson and M. Fabrizio, for useful discussions. R. 
Assaraf collaborated with this project in its earliest stage. This work was partly 
supported by the European Union, contract ERBFMRXCT970155 (TMR Fulprop), 
and by MIUR COFINOl, MURST COFIN99. 



7 



[1] V. Brouet, H. AUoul, F. Quere, G. Baumgartner, and L. Forro, Phys. Rev. Lett. 82, 
2131 (1999). 

[2] A. P. Ramirez, Supercond. Review 1, 1 (1994). 
[3] O. Gunnarsson, Rev. Mod. Phys. 69, 575 (1997). 

[4] J. H. Schon, Ch. Kloc, R. C. Haddon, and B. Batlogg, Science 288, 656 (2000). 
[5] J. H. Schon, Ch. Kloc, and B. Batlogg, Nature 408, 549 (2000). 
[6] M. Fabrizio and E. Tosatti, Phys. Rev. B 55, 13465 (1997). 

[7] M. Capone, M. Fabrizio, P. Giannozzi, and E. Tosatti, Phys. Rev. B 62, 7619 (2000). 
[8] P. J. Benning, F. Stepniak, and J. H. Weaver, Phys. Rev. B 48, 9086 (1993). 
[9] V. Brouet, H. AUoul, T. N. Le, S. Garaj, L. Forro, Phys. Rev. Lett. 86 4680 (2001). 
[10] M. J. Rosseinsky, D. W. Murphy, R. M. Fleming, and O. Zhou, Nature 364, 425 (1993). 
[11] T. Takenobu, T. Muro, Y. Iwasa, and T. Mitani, Phys. Rev. Lett. 85, 381 (2000). 
[12] T. T. Obu, H. Shimoda, Y. Iwasa, T. Mitani, M. Kosaka, K. U. Tanigaki, C. M. Brown, 

and K. Prassides, Mol. Cryst. Liq. Cryst. 340, 599 (2000). 
[13] H. Tou, N. Muroga, Y. Maniwa, H. Shimoda, Y. Iwasa, and T. Mitani, Physica B 281, 
1018 (2000). 

[14] O. Zhou, T. T. M. Palstra, Y. Iwasa, R. M. Fleming, A. F. Hebard, P. E. Sulewski, D. 

W. Murphy, and B. R. Zegarski, Phys. Rev. B 52, 483 (1995). 
[15] S. Margadonna, K. Prassides, H. Simoda, Y. Iwasa, and M. Mezouar, Europhys. Lett. 

56, 61 (2001). 

[16] O. Gunnarsson, E. Koch, and R. M. Martin, Phys. Rev. B 54, R11026 (1996). 

[17] E. Koch, O. Gunnarsson, and R. M. Martin, Phys. Rev. B 60, 15714 (1999). 

[18] The origin of the band splitting at the F point is twofold; (i) the hopping between 
neighbor fuUerene balls is different for pairs of molecules separated by different distances, 
along different crystal directions; (ii) different orbitals in the molecular tiu manifold 
composing the conduction band acquire a different energy according to their orientation 
relative to the neighboring Cg^, NH3 and K"*" species, forming the orthorhombic local 
environment. 

[19] A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. Mod. Phys. 68, 13 
(1996). 

[20] A. I. Lichtenstein, M. I. Katsnelson, and G. Kothar, Phys. Rev. Lett. 87, 067205 (2001). 

[21] M. Capone, M. Fabrizio, E. Tosatti, Phys. Rev. Lett. 86, 5361 (2001). 

[22] N. Manini, E. Tosatti, and A. Auerbach, Phys. Rev. B 49, 13008 (1994). 

[23] I. M. Lifshitz, Sov. Phys. JETP 11, 1130 (1960). 

[24] M. I. Katsnelson and A. V. Trefilov, Phys. Rev. B 61, 1643 (2000). 

[25] M. J. Rozenberg, Phys. Rev. B 55, R4855 (1997). 

[26] M. J. Rozenberg, G. Kotliar, and X. Y. Zhang, Phys. Rev. B 49, 10181-10193 (1994). 
[27] M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 (1994). 
[28] R. Bulla, Phys. Rev. Lett. 83, 136 (1999). 

[29] D. P. Arovas and A. Auerbach Phys. Rev. B 52, 10114-10121 (1995). 
[30] B. Sutherland, Phys. Rev. B 12, 3795 (1975). 

[31] Y. Q. Li, Michael Ma, D. N. Shi, and F. C. Zhang, Phys. Rev. Lett. 81, 3527 (1998). 

[32] Our implementation of the DMFT follows closely that of Ref. [Q. The chemical po- 
tential is adjusted at each iteration, so that the desired number of electrons per site 
is obtained at convergence We find advantageous to use a logarithmic mesh for 
the imaginary frequencies, with a low-frequency cutoff of the order of 0.02 W, to insure 
convergence with the rather small number of baths per band we can include (rib ~ 4-^5). 
Our results are in substantial accord with those of the traditional linear mesh. 

[33] Working at zero temperature is a quite important feature in our specific case, where 
the competition between the splitting A and the small quasi-particle residue z might 
be blurred, within a standard Monte Carlo approach, by finite temperature effects. The 
exact diagonalization method suffers, however, from a severe constraint on the number 
of discretized conduction band states which can be included. As a consequence, for 
realistic calculations with three or more degenerate orbitals a Monte Carlo approach is 
more appropriate [ p^ . 

[34] The package, togheter with the ultrasoft pseudopotent ials of C, N, H and local potential 

of K, can be found at the URL http: / /www.pwscf .or 
[35] S. Margadonna, K. Prassides, H. Shimoda, T. Takenobu, Y. Iwasa, Phys. Rev. B 64 



8 



132414 (2001). 

[36] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1974). 
[37] V. P. Antropov, O. Gunnarsson, and O. Jepsen, Phys. Rev. B 46, 13647 (1992). 
[38] R. L. Martin and J. P. Ritchie, Phys. Rev. B 48, 4845 (1993). 
[39] E. Koch, O. Gunnarsson, and R. M. Martin, Phys. Rev. Lett. 83, 620 (1999). 
[40] A. Goldoni, L. Sangaletti, F. Parmigiani, G. Comelli and G. Paolucci, Phys. Rev. Lett. 
87, 076401 (2001). 

[41] Polarization screening of the on-site interaction U is different in the two compounds, 
with poorer screening in NH3K3C60 leading to a larger effective U/W. However, a crude 
estimate suggests that the reduced screening of Ceo due to the slightly larger cell volume 
is probably over-compensated by the extra screening due to ammonia. As a net result, 
we expect the effective U to be substantially the same in the two compounds, as indeed 
seems to be the case in different fullerides. 

[42] H. Kitano, R. Matsuo, K. Miwa, A. Maeda, T. Takenobu, Y. Iwasa, and T. Mitani, 
Phys. Re v. Lett. 88, 096401 (2002). 

[43] E. Koch, |cond-mat/0 112329 . 

[44] S. Wehrh, D. Poilblanc, and T. M. Rice, Eur. Phys. J. B 23, 345 (2001). 
[45] S. Wehrli, private communication. 

[46] T. Costi and N. Manini, J. Low Temp. Phys. 126, 835 (2002). 



9 



2-band insulator 




AW 

FIG. 1. Qualitative zero-temperature phase diagram for the two-band Hubbard model at 
quarter filling (one electron per site) in the (7- A plane, where A is the anisotropy splitting 
of the two orbitals. The various phases and lines are described in Sect. The multicritical 
points M and M' are not necessarily distinct. 
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FIG. 2. DMFT zero-temperature phase diagram for the two-band Hubbard model at 
quarter fiUing (one electron per site), obtained by the exact diagonalization method in the 
paramagnetic sector. The discretization of the conduction band uses here nt = 4 bath states 
per orbital degree of freedom. 
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co/W 

FIG. 3. Spectral density Aa{uj) = ~Tr~'^lmGa{uj) at U/W = 1.5 across the Mott-Hub- 
bard transition for increasing anisotropy splitting A. ui is referred to the Fermi energy. Solid 
lines refer to the minority orbital a = 2. The multi-peak structures of the high-energy side 
bands are artifacts of the finite discretization of the conduction band. 
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All visible atoms of the central unit cell are shown together with some atoms in the neighbor 
unit cells. In NH3K3C60 the c-axis is in the 2 direction. 
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FIG. 5. Band structure of NHaKaCeo (a), obtained by an DFT-LDA calculation, com- 
pared to that of KsCfio (c). The corresponding densities of states are given in panels (b) 
and (d). In panel (b) the density of states of KsCeo (dashed line) is re-drawn for direct 
comparison with that of NH3K3C60 (solid line). 
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MOTT-HUBBARD INSULATORS 
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FIG. 6. A schematic (17, A) phase diagram for the d = 3-bands family of compounds 
KsCeo and NHsKsCeo- The dashed and dot-dashed arrows indicate two plausible paths of 
pressure-induced metallization of NHsKsCso- 
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